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ABSTRACT 

We present near- and mid-infrared photometry and spectroscopy from PAIRI- 
TEL, IRTF, and Spitzer of a metallicity-unbiased sample of 117 cool, hydrogen- 
atmosphere white dwarfs from the Palomar-Green survey and find five with excess 
radiation in the infrared, translating to a A.3tf 7 2 % frequency of debris disks. This 
is slightly higher than, but consistent with the results of previous surveys. Using 
an initial-final mass relation, we apply this result to the progenitor stars of our 
sample and conclude that 1 — 7M stars have at least a 4.3% chance of hosting 
planets; an indirect probe of the intermediate-mass regime eluding conventional 
exoplanetary detection methods. Alternatively, we interpret this result as a limit 
on accretion timescales as a fraction of white dwarf cooling ages; white dwarfs 
accrete debris from several generations of disks for ~10Myr. The average total 
mass accreted by these stars ranges from that of 200km asteroids to Ceres-sized 
objects, indicating that white dwarfs accrete moons and dwarf planets as well as 
Solar System asteroid analogues. 

Subject headings: infrared: planetary systems — infrared: stars — white dwarfs 

1. INTRODUCTION 

The study of extrasolar planetary systems has become a very active field that will 
continue to grow in the coming years and decades. Over the last two decades, several 
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hundred exoplanets have been discovered, proving planet formation to be a robust process 
in the Galaxy. The vast majority of known exoplanets have been discovered using just 
two techniques; namely, radial velocity and transits. To date, radial velocity studies have 
revealed the presence of well over five hundred planets and the number of transiting systems 



has passed the two hundred mark (see http : / / exoplanet.eu] and Schneider et al. 2011). 

A cursory examination of these planet-hosting stellar systems reveals a curious trend. 
The frequency of Jovian planets increases from 3% for M dwarfs to 14% for 2M Q stars 



(Johnson et al. 2010 Bowler et al. 2010). In addition, the distribution of known exoplanets 
as a function of host stellar mass has a strong cutoff at 3M , above which no exoplanets have 
been found. Could it be that planets do not form around intermediate-mass stars, or is this 



a selection effect? The planet formation models of Kennedy & Kenyon (2008) suggest the 
latter. They predict that the fraction of stars with giant planets shows a steady increase with 
mass up to 3M . In addition, the mass of the planets and the width of the regions where 
they form are predicted to increase with stellar mass. Probing intermediate-mass stellar 
systems for planets could help test these planet formation models. Since most exoplanet 
detection methods tend to break down in the intermediate-mass range, one can look at 
these systems in post main-sequence where the contrast between the planet's photometric 
signal and that of its host is greatly decreased. This work examines a sample of white dwarfs 
decedent from 1 — 7M stars and constrains the frequency of planetary systems in the elusive 
intermediate-mass regime. 



If planets survive post-MS evolution, they should be detectable around WDs (Burleigh 



et al. 2002; Gould & Kilic 2008). Unfortunately, there are still no confirmed planets around 



single WDs. The candidate planet around the pulsating WD GD 66 (Mullally et al. 2007a) 
is currently disputed (J. J. Hermes 2012, private communication). Perhaps an easier way to 
detect remnant planetary systems around WDs is to look for the tidally disrupted remains 
of exoplanets and moons in the form of circumstellar debris disks ( Debes fc Sigurdsson||2002 



Jura 


2003 


Kilic et al. 


2006) 



There are now two dozen WDs known to host dust and/or 
gas disks. Photospheric abundance analyses of these WDs show that the accreted metals 
originate from tidally disrupted minor bodies similar in composition to that of bulk Earth 



(Zuckerman et al. 2007 Klein et al. 2010, 2011 Dufour et al. 2010, 2012). Since at least 



one giant planet is required to perturb minor bodies out of their stable orbits (Debes et al. 



2012), photospheric pollution as well as circumstellar debris disks serve as tracers for giant 
planets at WDs. Therefore, the frequency of debris disks around WDs can be taken as a 
lower limit on the frequency of planets around WD remnants of intermediate-mass stars. 



The search for debris disks around WDs begins as a search for excess emission in the in- 
frared. Over the last two decades, several hundreds of WDs have been surveyed for infrared 
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excess due to substellar companions or circumstellar disks. The first dusty WD, G29-38, 



was identified by Zuckerman & Becklin (1987), see also Graham et al. (1990). Zuckerman 



& Becklin (1992) examined 200 WDs for excesses in the K-band, while Farihi et al. (2005) 



surveyed 371 WDs of all spectral types and constrained the frequency of substellar compan- 



ions to less than 0.5%. The discovery of the second dusty WD, GD 362, (Kilic et al. 2005 



Becklin et al. 2005) lead to the realization that dusty WDs may be found more readily by 



looking at metal-rich WDs. 

Previous surveys for debris disks around WDs focused on metal-rich WDs since so far 
all known WDs with disks also show metal absorption lines in their optical and ultraviolet 
spectra; DAZ and DBZ spectral types. The Poynting-Robertson drag timescale is signifi- 



cantly shorter than the evolutionary timescales for WDs (von Hippel et al. 2007 Rafikov 



2011b), hence we expect every WD with dust /gas disks to accrete metals from these disks 
and appear metal-rich. Zuckerman et al. ( 2003[ 2010) find that about 30% of DA and DB 
WDs show trace amounts of metals in their photospheres. However, only ~ 20 — 30% of 



DAZs, the ones with the highest accretion rates, host dust (Kilic et al. 2008 Farihi et al. 



2009). The source of metals for the other DAZs remains uknown, but it is likely that the 



majority of DAZs accrete many small asteroids that do not form a large debris disk (Jura 



2008). 



Mullally et al. (2007a) survey 124 WDs with T eB = 5,000 - 170,000K using mid- 



infrared Spitzer photometry. Their search for infrared excess finds 2 dusty WDs for a disk 
frequency of 1.6%. However, this survey extends beyond the temperature range within which 
solid dust orbiting interior to the tidal radius can persist. Based on Spitzer observations of 



mostly metal-rich WDs, Farihi et al. (2009) estimate a debris disk frequency of 1% to 3% 



for WDs younger than 500Myr. Debes et al. (2011) recently searched for infrared excesses 



amongst an unbiased (in terms of metallicity) WD sample using the Wide-field Infrared 



Survey Explorer (WISE, |Wright et aL 2010) data. They found about 1 - 5% of these WDs 



to be dusty. However, the large beam size of WISE allows for many false positives due to 
background contamination, and many of these dusty WD candidates need follow-up higher 
spatial resolution infrared imaging data to confirm the observed infrared excess. 

Here we present a near- and mid-infrared photometric and spectroscopic survey of an 
unbiased (in terms of metallicity) sample of DA WDs from the Palomar-Green (PG) survey. 
Our survey does not focus on metal-rich WDs and thus constrains the disk frequency repre- 
sentative of the general WD population. We target 117 PG WDs with accurate temperature, 
mass, and age estimates (Liebert et al. 2005) and T e g = 9, 500 — 22, 500K. With 117 objects 



we sample temperatures at which solid dust can persist within the tidal radius and thus 



improve upon the constraint obtained by Mullally et al. (2007a) whose sample includes 69 
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WDs in this temperature range. The WDs in our sample are the descendants of MS stars 
with masses between 1 — 7M , so this study extends the search for planets well beyond the 
mass range available to conventional exoplanet detection methods. Section 2 describes our 
near- and mid-infrared photometry and spectroscopy data, while Section 3 presents the spec- 
tral energy distributions (SEDs) and the dusty WDs in the PG survey. Section 4 presents a 
discussion of the frequency of disks and planets around WDs and their progenitor MS stars. 



2. OBSERVATIONS 



2.1. The Sample 

Almost all known dusty WDs are relatively young with cooling ages of < lGyr and tem- 
peratures in the range 9, 500 — 22, 500K (with the exception of G166-58, Farihi et al. 2008). 
This is probably because a second heavy bombardment phase occurs around intermediate- 
mass stars in their post-MS evolution due to planetary migration around the mass-losing star 



(Debes et al. 2012 Bonsor & Wyatt 2012). We would expect to see a peak in the frequency 



of collisions, tidal disruptions, and disks around the younger WD remnants of such stars 



(Debes fe Sigurdsson 2002). 



Consider the disk around G29-38 as a typical system. Jura (2003) models the disk 
emission and obtains a ring of dust extending from R in = 0.1R Q to R out = 0.4 — O.9i? 
and T eff = 300 — 600K. Such a disk around a star twice or three times as hot as G29-38 
would be sublimated and invisible in the IR. Thats why our target selection using effective 
temperature gives a better measurement of the real frequency of disks. 

We select 117 apparently single DA WDs from the PG survey with T e g- = 9,500 — 
22, 500K, where we are most efficient in finding the disks using near- and mid-infrared data. 



Liebert et al. (2005) performed a detailed model atmosphere analysis of all DA WDs in 



the PG survey and provided temperature, surface gravity, mass, and age estimates. Their 
spectroscopy did not have enough resolution to detect the metal-rich DAZs among these 
targets; our survey does not have a metallicity bias, and hence it has the potential to reveal 
the true frequency of dusty disks without knowing the fraction of metal-rich WDs among 
the nearby WD population. 
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2.2. Near-infrared Photometry 

We obtained simultaneous JHK S imaging of 78 WDs using the Peters Automated Imag- 
ing Telescope (PAIRITEL, [Bloom et al^|200~6| ) between 2009 Feb and 2010 Apr. PAIRITEL 



is the old Two Micron All Sky Survey (2MASS, Cutri et al. 2003) telescope and it uses the 
same camera and filter set. The total exposure time for our targets ranged from 70 s to 
1,500 s depending on the expected brightness of each target. We use the PAIRITEL pipeline 
version 3.3 processed images and standard IRAF DAOPHOT routines to perform aperture 
photometry on every 2MASS source detected in the images. We use the 2MASS stars to 
calibrate the photometry for each object. 

Table 1 presents the PAIRITEL photometry (and the 2MASS photometry of a few 
relatively bright targets) for our sample along with the physical parameters obtained from the 



model atmosphere analysis of Liebert et al. (2005). We compare our PAIRITEL photometry 



to that of 2MASS in Figure [T} Only about half of our targets are detected by 2MASS in 
the K-band, and usually with large photometric errors. Our survey goes significantly deeper 
than 2MASS, and provides improved photometry with smaller errors, which is essential 
for identifying objects with slight K— band excesses. This figure demonstrates that the 
PAIRITEL photometry agrees fairly well with the 2MASS photometry and that it can be 
used reliably to constrain the near-infrared SEDs of our targets. 

Figure [2] shows J — H and J — K colors of our sample of stars as a function of their 
effective temperatures. The colors for the majority of our targets follow the model predictions 
of J — H = J — K m mag, though with some scatter. There are four WDs with J — H 
excesses, though with relatively noisy 2MASS photometry. However, from the available mid- 
infrared photometry from Spitzer and/or WISE, we conclude that there is no evidence of 
excess in these objects. The J — H photometry for the rest of the targets are consistent 
with the emission from single WDs. Several of these apparently single WDs have positive 
J — K colors, suggesting extra emission from a cool companion or a circumstellar debris 
disk. We chose to further examine those objects with positive J — K colors by obtaining 
near-infrared spectra and mid-infrared photometry. Objects with apparently red colors were 
targeted however, depending on the timing of the observing runs and conditions, some of 
the objects with blue colors were also followed up as a sanity check. 



2.3. Near-infrared Spectroscopy 



We obtained low-resolution near- infrared spectra of 41 WDs over several nights in 2011 
April, August, and September using the 3-m NASA InfraRed Telescope Facility (IRTF) 
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equipped with the 0.8 — 5.4 Micron Medium-Resolution Spectrograph and Imager (Spex; 



Rayner et al. 2003). The observing setup and procedures are similar to those of Kilic et al. 



(2012). We use a 0.5" slit to obtain a resolving power of 90-210 over the 0.7 — 2.5/xm range. 



The observations are taken in two different positions on the slit separated by 10". We use 
internal calibration lamps (a 0.1W incandescent lamp and an Argon lamp) for flat-fielding 
and wavelength calibration, respectively. To correct for telluric features and flux calibrate 
the spectra, we use the observations of nearby A0V stars at a similar airmass to the target 



observations. We use the IDL-based package SPEXTOOL version 3.4 (Cushing et al. 2004) 
to reduce the data. 

To demonstrate that our near-infrared observations and reductions are reliable, we also 
obtained a near-infrared spectrum of the previously known dusty WD GALEX J193156. 8+01 1745 
(hereafter J1931+0117; |Vennes et aLl|20lo| |Debes et al-lpoTTj |Melis et al.|[20TT| ). Figure g 
shows the near- infrared photometry and our IRTF spectrum of this object. There are two 



nearby red sources within 2" of the WD (Melis et al. 2011). Our IRTF observations were 



obtained under average seeing conditions of 0.8" and our 0.5" slit minimizes the contamina- 
tion from the nearby sources. A slight excess is clearly detected redward of 2.1/im, typical of 



dusty WDs (Kilic et al. 2006). Melis et al. (2011) also obtained a near-infrared spectrum of 



J1931+0117 on the Magellan 6.5m Baade telescope, which shows a significantly larger excess 
redward of 1.6/mi. Those authors discuss the potential problems with the flux calibration 
of their spectra. They model the excess around J1931+0117 with a flat disk model with 
inner and outer temperatures of 1400 and 1200 K, essentially a narrow ring of dust. Our 
IRTF spectrum demonstrates that the excess around J1931+0117 mostly shows up redward 
of 2.1/zm and it is similar to the other dusty WDs with extended (up to IRq) disks. Of 
course, the outer edge of the disk around J1931+0117 is currently unconstrained due to 
the lack of uncontaminated mid-infrared photometry. Regardless of these issues, this figure 
demonstrates that we are able to identify dusty WDs even with slight K— band excesses. 



2.4. Spitzer Photometry 



We used the warm Spitzer equipped with the InfraRed Array Camera (IRAC; Fazio 



et al. 2004) to obtain infrared photometry of 11 WDs between 2010 Sep and 2011 Feb for 



program number 70023. Based on our preliminary analysis of the PAIRITEL data, these 11 
sources appeared to have slight K— band excesses indicative of debris disks. We obtained 3.6 
and 4.5/xm images with integration times of 30 s or 100 s for nine dither positions. We use 
the IDL ASTROLIB packages to perform aperture photometry on the individual corrected 
basic calibrated data frames from the S18.18.0 pipeline reduction. 
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Following the IRAC calibration procedures, we correct for the location of the source in 
the array before averaging the fluxes of each of the dithered frames at each wavelength. We 
also correct the Channel 1 (3.6/xm) photometry for the pixel phase dependence. We estimate 
the photometry error bars from the observed scatter in the nine images corresponding to the 



dither positions. We add the 3% absolute calibration error in quadrature (Reach et al. 



2005). Finally, we divide the estimated fluxes by the color corrections for a Rayleigh- Jeans 



spectrum, except for the dusty WDs in our sample. Table [2] presents our Spitzer IRAC 
photometry for 11 PG WDs. 



3. RESULTS 

Our near-infrared photometric observations revealed several sources with potential K — band 
excesses. We followed up 11 of the most interesting targets with Spitzer IRAC and 41 stars 
in total at the IRTF. Figure [4] shows the SEDs of the 41 stars observed at the IRTF. Here, 
and in the rest of the figures, the PAIRITEL and 2MASS photometry are shown as blue and 
green points, respectively. The black lines show the observed spectra and the red solid lines 
show the predicted photospheric emission for each star assuming a blackbody. The IRTF 
spectra display telluric correction problems around 1.4 and 1.9/zm, but otherwise they closely 
follow the predicted blackbody distributions for the majority of the targets in our sample. 
PG 0048+202 and PG 1720+361 are two of the stars for which our PAIRITEL photometry 
indicated K— band excesses, but our IRTF spectra show that there are no significant ex- 
cesses in the K— band for these two stars. There is only one WD in our IRTF sample where 
a significant K— band excess is detected, PG 1541+651. The observed flux stays essentially 
constant between 2.0 and 2.5 /im, which is typical of dusty WDs. The detection of Ca II 3933 
A absorption in the HIRES spectrum of this object (Xu et al. in preparation) corroborates 
the presence of circumstellar debris. 

Figure [5] shows the near- and mid- infrared data on 11 WDs that we targeted with 
Spitzer, including PG 1541+651. Eight of these stars also have IRTF spectroscopy available. 
Our Spitzer observations confirm the results from the IRTF observations; only PG 1541+651 
shows a significant infrared excess that is compatible with a debris disk (see |Kilic et al.|20l2| . 



The SEDs for the other ten objects in this figure are consistent with photospheric emission 
from each star. 

There are four other previously known dusty WDs in our sample; PG 1015+161, PG 
1116+026, PG 1457-086, and PG 2326+049. Figure © displays the SEDs for these four 
WDs plus PG 1541+651. PG 2326+049 has the most prominent excess, while the excess at 
PG 1457—086 is the most subtle. The variation in the strength of these excesses is likely 
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due to a variation in the disk geometry from object to object. There is another known dusty 
WD in the |Liebert etaL] ( pOOo) DA WD sample, PG 1456+298 (G 1 66-58, [Farihi et al.|[2008| ) 



that is cooler than the lower temperature limit of our sample and hence it is not included in 
Figure |6j 

Out of the 117 DA WDs in our sample, the majority have SEDs that closely follow a 
single blackbody curve, while five WDs show a significant deviation from their blackbody 
model in the near- to mid-infrared (see Figure [6]). In light of the recent WISE All-Sky 
Data Release, we examine our sample of objects for infrared excesses in the WISE data 
and plot these results in Figure [7} Four of the five dusty WDs in our sample clearly show 
infrared excesses in the WISE bands. The fifth dusty WD, PG 1457—086, is also detected 
in the WISE observations, but with large photometric errors. There is another target, PG 
1519+384, with positive H — Wl and H — W2 colors. However, there are two nearby sources 
within several arcseconds of this target in the Sloan Digital Sky Survey images. The WISE 
photometry for PG 1519+384 is likely contaminated by these background sources. Hence, 
the WISE data do not reveal any new dusty WDs in our sample other than the five systems 
known. 



Farihi et al. (2005) finds a frequency of less than 0.5% for brown dwarf companions 
of WDs. Given our sample size of 117 we would expect to find less than one brown dwarf 
companion in our sample. It is not surprising, then, that we do not find any substellar 
companions. Through the initial selection process, we excluded objects with potential M- 
dwarf companions. 



DISCUSSION 



4.1. The Frequency of Debris Disks Around WDs 



Out of the 348 DA WDs analyzed by Liebert et al. (2005), 308 are apparently single. Six 



of these PG WDs are known to host dust disks, corresponding to a debris disk frequency of 
1.9%. This is similar to the frequency of disks, 1.6%, around the nearby bright WD sample 



of Mullally et al. (2007a). Of course, these estimates ignore the fact that most disks occur 



around young WDs (see |Kilic et al.||2009bp . The process that creates the debris disks around 
WDs seem to be more efficient at younger ages. 

We use the binomial probability distribution to compute the upper and lower limits on 
the frequency (p) of disks. The binomial distribution function gives the discrete probability 
distribution of obtaining exactly n successes out of N trials (where the result of each trial 
is true with probability p and false with probability 1 — p). The probability, P n (N,p), that 
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a survey of A" WDs will detect n debris disks, when the true frequency of disks is p is given 



by (Burgasser et al. 2003 Appendix 



Pn(N, P ) 



N 



n 



p n {l-p) 



N-n 



N\ 



n\(N - n)\ 



p n (l-p) 



N-n 



X 



Since this probability function is not symmetric about its maximum value, we report the 
range in p that delimits 68% of the integrated probability function, equivalent to la Gaussian 
limits. For A" = 117 and n — 5, there is a 68% chance of p being between 3.1% and 7% 
(see Figure |8j). Thus, p = 4.3tfl%- This frequency estimate is slightly higher than, but 



consistent with, the 3% disk frequency derived by Farihi et al. (2009). 



4.2. Planets Around Intermediate-Mass Stars 



The 117 WDs in our sample have an average mass of O.6M and mass range of 0.4 — 
1.2M . We estimate the progenitor MS masses using the initial-final mass relation derived 



by Kalirai et al. (2008) and Williams et al. (2009). We obtain an average progenitor mass 



of 2.4M Q and mass range of 1.1 — 7.2M (see Table [T]) where we neglect WDs below O.5M . 
Based on the observed frequency of disks around our WD targets, we find that the probability 
of finding planetary systems around their progenitor 1 — 7M MS stars is at least 4.3^^2%. 
This result serves as an indirect probe of the intermediate-mass stellar systems which have 



thus far challenged conventional planetary detection techniques (Gould & Kilic 2008). 



The dusty WDs in our sample, PG 1015+161, 1116+026, 1457-086, 1541+651, and 
2326+049 are the descendants of 2.1 — 2.8M MS stars. At first glance, this argues that 
planet formation may be most efficient around 2 — 3M stars. However, given the age and 
star formation history of the Galactic disk, our sample is dominated by the ~ O.6M WD 
remnants of 2 — 3M MS stars. In fact, 80% of our sample consists of WD descendants 
of 1 — 3M MS stars. Hence, the lack of discoveries of dusty disks, and therefore remnant 
planetary systems, around the massive WD remnants of M > 3M stars may be due to 
small number statistics. There are 22 systems in our sample with M > 3M progenitors, 
we would expect to find 1 dusty WD in a sample of 22 stars. The probability of finding 
zero dusty WDs in a sample of 22, when the expected number is one, is 30% (see Equation 
1). Hence, larger surveys of massive WDs are required to detect disks around them and to 
constrain the frequency of planets around their intermediate-mass progenitor MS stars. 
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4.3. Disk Lifetimes 



An alternative interpretation of the 4.3% disk frequency around WDs is that perhaps 
100% of WDs have remnant planetary systems that can create debris disks, but the disks only 
persist for 4.3% of the cooling age of the WD. For the average cooling age of our sample of 
WDs, 300Myr, this implies that disks survive around typical WDs for ~10Myr. We note that 
this is the total time a WD accretes from circumstellar debris disks, even though individual 



disks may last for a significantly shorter amount of time. Kilic et al. (2008) and Farihi et al. 
(2009) estimate typical disk lifetimes of ~ 10 5 yr. Hence, up to 100 tidal disruption events 



that create circumstellar debris disks may occur for any given WD. 



von Hippel et al. (2007) use a geometrically thin, optically thick disk model to determine 
typical accretion rates of WDs. Their result (~ 10 9 g s _1 ) means a WD with the disk lifetime 
typical of our sample will accrete a total mass equivalent to a 400 km Solar System asteroid 
(assuming an asteroid density of 3 g cm -3 ). Rafikov (2011a) demonstrate that the Poynting- 
Robertson drag can explain accretion rates of around 10 8 g s _1 from circumstellar debris 
disks onto the WDs. Taken at face value, this implies a total accreted mass of 10 22 g, a 
200km asteroid, over the disk-hosting lifetime of a WD. However, dusty WDs show accretion 



rates of up to 10 11 g s 1 (Dufour et al. 2012). To explain the observed high accretion rates, 



Rafikov (2011b) and Metzger et al. (2012) propose a runaway accretion scenario where rapid 



transport of metals from the disk results from interaction between spatially coexisting dust 
and gas disks. In this scenario, a WD accretes 10 22 g of metals in 10 5 years. Hence, typical 
WDs may accrete up to 10 24 g of metal over lOMyr, the total time they host circumstellar 
disks. This is comparable to the mass of the dwarf planet Ceres and Pluto's moon Charon. 
Of course, if the tidally disrupted object is icy, the total accreted mass may be significantly 



higher (although see Jura & Xu 2012). Therefore, the debris surrounding WDs is likely 
supplied through disruption of Solar System asteroid analogues, moons, and dwarf planets. 



5. CONCLUSIONS 

We present near- and mid-infared observations and a comprehensive study of the SEDs 
of a metallicity- unbiased sample of 117 DA WDs from the PG survey. Only five of our targets 
show excess radiation from debris disks, indicating a debris disk frequency of 4.3+ 2 ;2% which 
fits in well with previous surveys of this kind. We interpret this frequency as a lower limit 
to the frequency of planets around the 1 — 7M progenitor MS stars; an indirect result 
for the intermediate-mass regime in which conventional exoplanetary detection methods are 
insensitive. Alternatively, we interpret the observed frequency of disks as the fraction of time 
a WD accretes from its debris disks over the entire cooling age of the star. We estimate that 
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typical WDs may accrete metals from several generations of disks for ~10Myr, corresponding 
to accretion of metals up to 10 24 g. This means WDs are capable of accreting bodies as large 
as dwarf planets as well as Solar System asteroid analogues. 
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Table 2. Spitzer IRAC Photometry 





3.6 /im 


4.5 /zm 


PG 




(^Jy) 


0826+455 


268.6 ± 8.7 


172.0 ± 5.9 


0956+021 


88.6 ± 3.5 


53.7 ± 2.8 


1036+086 


42.5 ± 1.5 


26.6 ± 1.3 


1122+546 


118.2 ± 4.3 


77.6 ± 3.1 


1149+058 


248.0 ± 8.1 


163.4 ± 5.7 


1325+168 


64.7 ± 2.2 


42.3 ± 1.7 


1541+651 


435.4 ± 13.5 


483.9 ± 15 


1605+684 


60.9 ± 2.7 


41.2 ± 2.2 


1720+361 


155.9 ± 5.4 


97.6 ± 3.7 


2303+243 


170.8 ± 5.9 


110.1 ± 4 


2306+131 


146.2 ± 5.2 


92.0 ± 3.5 



Note. - Spitzer Cycle 7 photometry 
fluxes in IRAC Channels 1 (3.6 /im) and 2 
(4.5 ttm). 
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Fig. 1.— The PAIRITEL vs. 2MASS photometry for our sample. The PAIRITEL data 
agrees fairly well with the 2MASS photometry. 
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Fig. 2. — J — H and J — K colors vs. temperature for our sample of DA WDs from the 
PG survey. Red squares and blue triangles show the confirmed dusty WDs and objects with 
follow-up IRTF spectroscopy, respectively. Solid lines show the predicted sequences for DA 
WDs (|Bergeron et~aL]|2011|>. 
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Fig. 3. — Spectral energy distribution of the known dusty WD GALEX J1931+0117. Green 



points represent 2MASS photometry, while blue points represent the photometry from Melis 



et al. 2011 The IRTF spectrum is shown in black and the red line represents the predicted 



contribution from a blackbody stellar photosphere. The noise near 1.9 and 2.5 /im falls in 
regions of poor telluric correction. The 1.28 fim feature is Pa/5 absorption. The detection 
of slight K-band excess around this target demonstrates that our near-infrared observations 
and reductions are reliable. 
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Fig. 4. — Spectral energy distributions of 41 WDs observed at the IRTF. The PAIRITEL 
and 2MASS photometry are shown as blue and green points, respectively. The black lines 
show the observed spectra and the red solid lines show the predicted photospheric emission 
for each star assuming a blackbody. The IRTF spectra display telluric correction problems 
around 1.4 and 1.9 /im. 



(p ) 

Fig. 5. — SEDS of the 11 WDs with new Spitzer photometry (shown in purple). The symbols 
are the same as in Figure [3j The IRTF spectra display telluric correction problems around 
1.4 and 1.9 fim. Our Spitzer observations confirm the results of the IRTF observations; only 
PG1541+651 shows a significant infrared excess. 
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Fig. 6. — SEDs of the five dusty WDs in our sample. The symbols are the same as in Figure 
[5| The IRTF spectra display telluric correction problems around 1.4 and 1.9 fim. 
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Fig. 7. — Color-color diagrams using PAIRITEL/2MASS and WISE photometry for our 
WD sample. Red triangles mark the dusty WDs. Four of the five dusty WDs in our sample 
clearly show infrared excesses in the WISE bands. The fifth dusty WD, PG 1457—086, is 
also detected in the WISE observations, but with large photometric errors. 
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Fig. 8. — Probability function for WDs with debris disks. Probability, P n (N,p), of finding 
n = 5 debris disks as a function of assumed true debris disk frequencies, p, where N = 117. 
Dashed lines delineate the region containing 68% probability, equivalent to la Gaussian 
limits. The solid line indicates the peak of the distribution. The derived frequency is 



